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ABSTRACT 

Context. In a recent study, star-forming galaxies with He n ,11640 emission at moderate redshifts between 2 and 4.6 have been found 
to occur in two modes that are distinguished by the width of their Hen emission lines. Broad Hen emission has been attributed to 
stellar emission from a population of evolved Wolf-Rayet (WR) stars. The origin of narrow He n emission is less clear but has been 
attributed to nebular emission excited by a population of very hot Pop III stars formed in pockets of pristine gas at moderate redshifts. 
Aims. We propose an alternative scenario for the origin of the narrow Hen emission, namely very massive stars (VMS) at low 
metallicity (Z), which form strong but slow WR-type stellar winds due to their proximity to the Eddington limit. 

Methods. We estimated the expected He n line fluxes and equivalent widths based on wind models for VMS and Starburst99 popula¬ 
tion synthesis models and compared the results with recent observations of star-forming galaxies at moderate redshifts. 

Results. The observed He n line strengths and equivalent widths are in line with what is expected for a population of VMS in one or 
more young super-clusters located within these galaxies. 

Conclusions. In our scenario the two observed modes of He n emission originate from massive stellar populations in distinct evo¬ 
lutionary stages at low Z (~ 0.01 Z 0 ). If this interpretation is correct, there is no need to postulate the existence of Pop III stars at 
moderate redshifts to explain the observed narrow He n emission. An interesting possibility is the existence of self-enriched VMS 
with similar WR-type spectra at extremely low Z. Stellar He n emission from such very early generations of VMS may be detectable 
in future studies of star-forming galaxies at high redshifts with the James Webb Space Telescope (JWST). The fact that the Hen 
emission of VMS is largely neglected in current population synthesis models will generally affect the interpretation of the integrated 
spectra of young stellar populations. 

Key words. Stars: Wolf-Rayet - stars: mass-loss - stars: Pop III - galaxies: starburst - galaxies: stellar content - galaxies: star 
clusters: general. 


1. Introduction 


A fundamental question in Astronomy is the question of the 
sources of the “First Light” ending the Cosmic Dark Ages and 
beginning the process of reionisation. The James Webb Space 
Telescope (JWST) promises direct access to this critical period 
via observations of the first star-forming galaxies at redshifts 
z > 10. Dual Lycr ,11216 and Hen ,11640 emission is believed 
to be the main indicator of the first (Pop III) stars formed from 
pristine gas in these objects (e.g. lJohnson et al]l2009 ). The rea¬ 
son is that only at extremely low metallicities (Z < 10 5 ), young 
massive stars are believed to be hot enough to excite He n in their 
surrounding Hu regions (e.g. !Schaereril2002l 120031) . The inves¬ 
tigation of star-forming galaxies with He n emission at moderate 
redshifts, which are accessible with current ground-based instru¬ 
mentation, is thus an important preparation for future studies of 
the first star-forming galaxies with the JWST. 


Such a study has recently been performed by ICassata et al.1 
( 2013 -). who identified different groups of He n emitters with nar¬ 
row and broad emission lines in a sample of star-forming galax¬ 
ies at redshifts 2 < z < 4.6. While the broad emission is usu¬ 
al ^attributed to stellar emission from Wolf-Rayet (WR) stars, 
ICassata et al.l interpreted the narrow emission as being due to 
ionising radiation from Pop III stars that must have formed in 


pockets of residual pristine gas at relatively low redshift, or a 
stellar population that is rare at z ~ 0 but more common at z ~ 3. 

In this letter we discuss the possibility that the narrow He n 
emission might also be of stellar origin. To this purpose we dis¬ 
cuss in Sect.[2]how the spectral properties of WR-type emission 
line stars are expected to change at low Z, estimate the expected 
line shapes and fluxes in Sect. [3] and compare our results with the 
observations of ICassata et al.l in Sect. [4] Conclusions are drawn 
in Sect.|5] 

2. Hen emission from WR stars at different 
metallicities Z 

The emission-line spectra of WR stars are formed in their strong, 
optically thick stellar winds. As a result of the large optical 
depth of these winds, ions like He hi recombine already within 
the wind, leading to recombination cascades with strong emis¬ 
sion lines between excited levels of the subordinate ions. Promi¬ 
nent examples are Hen41640 and Civ41550 in the UV, and 
He ii 44686 and C iv 45808 in the optical wavelength range. The 
width of the emission lines is mainly determined by the terminal 
wind speed of the stellar wind. To understand why we natu¬ 
rally expect two populations of WR stars with broad and narrow 
emission lines at low Z, we need to discuss the origin of the dif¬ 
ferent types of WR-type stellar winds. 
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2.1. WR populations in the Galaxy and LMC 


Comprehensive studies of the local population of WR stars of 
the nitroge n sequence (W N stars) have been performed for the 
Galaxy by Hamannetah (2006) and fo r the Large Magellanic 
Cloud (LMC) by Hainich et al. (120141) . Both works identified 
two different groups of WN stars. The first group are classical 
WN stars with H-poor/H-free surface compositions and hot tem¬ 
peratures up to 140 kK, in agreement with a core He-burning 
evolutionary stage. The second group are the so-called WNh 
stars with relatively H-rich surface compositions and cooler 
temperatures, in agreement with a core H-burning evolutionary 
stage. The WNh stars are further characterised by remarkably 
high luminosities in excess of ~ 1O 59 L 0 , while the group of clas¬ 
sical WN stars has luminosities below this value. 

The dividing luminosity at roughly 1O 5,9 L 0 suggests that the 
classical WN stars are core He-burning objects in a post-red su¬ 
pergiant (post-RSG) phase. This picture is further supported by 
the work of lSander et alJ (120121) . who found that the WR stars of 
the carbon sequence (WC stars), which are showing the products 
of He-burning at their surface, form a natural succession of the 
sequence of classical WN stars in the HR diagram. 

In this work we are mainly interested in the group of WNh 
stars with luminosit i es in the range 1O 59 L 0 < L < 1O fl9 L 0 (cf . 


Hamannetal. 2006; Crowth er et alJ 12010: Hain ich et alJl2014t 
Bestenlehner et alJl2014l) . The HRD positions of these objects 


suggest that they are very massive stars (VMS) with initial 
masses in exce ss of ~ 100 M 0 , and reachi ng up to 300 M 0 (cf . 

I Vink et alll2015l) . In their LMC sample, lHainich et alJ (120141) 
found 12% of all putatively single WN stars in the WNh stage. 
Combined with their high luminosities, this suggests that WNh 
stars have substantial impact on the integrated He n emission of 
star-forming galaxies. 


2.2. WR mass loss at low metallicities Z 


Observational evidence for a Z-dependence of WR mass-loss 
has been found by comparisons of the sample of WC stars 
with k now n distances in the Galax y and LMC by IGrafener et alJ 
(Il998li and ICrowther et alJ ((2002). In the following, we invoke 
the results of theoretical studies to understand how the proper¬ 
ties of WR stars are expected to change over a broad range of 
Z. 


The winds of class i cal W R stars have been modelled 
by GiMengr_&Hamaim| (120051) (for early WC subtypes) and 
IVink & de Koterl (120051) (for late WC and WN subtypes). Al¬ 
though there are still problems to reproduce the winds of in - 
termediate spectral subtypes (cf. IGrafener & HamannI 120051) . 
these works suggest that the winds of WR stars are radiatively 
driven, predominantly by millions of spectral lines of the iron- 
group elements. As a consequence, the mass-loss rates of WR 
stars are expected to depend on the environment metallicity Z 
dVink & de Koteiil2005i) . 


The winds of luminous W Nh stars have been modelled by 
IGrafener & Hamaiml (120081) and lVink et alJ (1201 ll) . These works 
suggest that the proximity to the Eddington limit plays a key role 
for the formation of WR-type stellar winds. This has been fur¬ 
ther suppo rted in systemat i c stud ies o f samples of VMS in youn g 
clusters by Grafener et alJ (120 1 ll) andlBestenlehner et alJ d2014l) . 
According to Graf ener & HamannI d2008l) . Eddington factors of 
the order of T e ~ 0.5 are required to initiate WR-type stellar 
winds at solar metallicity Z 0 (here T e = (x e L)/(4ncGM) de¬ 
notes the classical Eddington factor that is defined with respect 
to the opacity of free electrons Xe)- The required Eddington fac¬ 


tors suggest that the masses of WNh stars are higher than those 
of core He-burning stars, meaning that WNh stars are indeed 
VMS in the phase of core H-burning. 

How does the spectral appearance of the two groups of WR 
stars change at low metallicities? As a result of their dependence 
on the iron-group opacities, the winds of classical WR stars are 
expected to be weaker at low Z dVink & de Koterll2005l). result- 
ing in significantly weaker emission lines dCrowther & Hadfieldl 
120061) . However, the situation is slightly different for WC stars 
because their surfaces are substantially self-enriched with the 
products of He-burning, namely C O, and most likely Ne and 
Mg. According to iVink & de Koterl d2005l) . these elements gain 
in importance for Z < 0.1 Z 0 , and may even dominate WR wind 
driving for Z < 10 3 Z 0 . WC stars may thus be the main source 
of broad He u emission from classical WR stars at low Z. Be¬ 
cause of the high carbon abundances of WC stars, this emission 
will most likely be accompanied by strong C iv emission. 

For the group of very massive WNh stars, the depen- 
de nce on the environment m etallicity Z has been discussed 
by Grafener & HamannI ( 2008 ). According to their models, the 
proximity to the Eddington limit is mandatory to support WR- 
type winds in this regime. The required Eddington factors in¬ 
crease for decreasing metallicity, from T e « 0.5 for Z = Z 0 , to 
r e ~ 0.85 for Z = O.O1Z 0 . As a consequence, the resulting ef¬ 
fective escape velocities decrease, and very low terminal wind 
velocities of only few lOOkm/s are predicted for WNh stars in 
this regime. At the same time, the mass-loss rates do not change 
significantly, resulting in significantly higher wind densities and 
much stronger and narrower WR emission lines than for classi¬ 
cal WR stars. Very massive WNh stars at low Z may thus be the 
sources of the narrow He n emission observed by ICassata et al.l 
( 20131) . In contrast to classical WC stars, they have very low car¬ 
bon abundances and are not expected to show detectable Civ 
emission. 


3. Expected Hen emission from populations of VMS 

Because of their short lifetimes (2-3 Myr), VMS are predom¬ 
inantly found in young massive clusters. To estimate the line 
fluxes expected from populations of VMS at cosmological dis¬ 
tances, we start with an examination of young massive clusters 
in our neighbourhood and their stellar He n emission in Sect. 13. II 
and extrapolate our results to low Z in Sect. 13. 21 

3. 1. VMS in young massive clusters 

IWoffordet al.l ( |2014|) recently found evidence for Hen emis¬ 
sion from VMS in the nearby starburst galaxy NGC3125 with 
roughly LMC metallicity. For the young cluster Al they had 
to invoke a population of VMS to explain the large observed 
Hen equivalent width (EW(Hen) ~ 7 A, cf. also IChandar et al] 
120041) as the cluster is too young (3-4 Myr) to host evolved WR 
stars. Based on different UV extinction corrections, they deter¬ 
mined an intrinsic He n line luminosity in the range log(LHen) = 

39.1.. .40.0 (in erg/s) and a cluster mass of log(M c i/M 0 ) = 

5.1.. .5.9. 

To investigate whether VMS can account for the large 
EW(Hen), we also investigated R136 in the centre of 30 Dor of 
the LMC, the most prominent starburst-like cluster in our neigh¬ 
bourhood that is spatially resolved. Using the integrated IUE 
spectrum of the central 5 pc from IVacca et ahl (t 19951) . we mea¬ 
sured a similarly high EW(He n) ~ 5A for this object. 
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For the sum of the seven known WNh stars in this region 
(five W N5h, oneWN6(h), and oneWN9ha; cf. ICrowther et alJ 
[20loh . we estimate in total log(LHen) ~ 38.6. This esti¬ 
mate is based o n a synthetic mode l for th e LMC WN5h star 
VFTS 682 from iBestenlehner et all (1201 ll) with log(L/L 0 ) = 
6.5 and log(LHen) = 37.5, which was scaled to the total lu¬ 
minosity of th e seven objects (log(L/L 0 ) — 7.58 + 0.0 3, as 
obtained from Crowther et alJ 1201 Ol lHainich et alJ 1201 4 and 
Best enlehner et al.ll20l4) . We note that this value is uncertain as 
for some of the stars different luminosities have been obtained in 
different works, and one of the stars, Mk 34, is a known binary. 

To estimate the intrinsic continuum flux, we used Star- 
b urst99 models dLeitherer et al. 1 17999) that are similar to those 
of lWofford et al.1 ( 2014 cf. their Sect. 3.5) but for Z = 0.008. For 
a cluster mass of log(M c i/M 0 ) ~ 4.7 dDoran et al.ll201 3if1 and an 
age of 2—4Myr, we obtain a continuum flux log(Luv) ~ 37.74 
(in erg s -1 A -1 ). Together with our estimate of log(Lneu) ~ 38.6, 
this results in EW(Hen) ~ 7 A, in good agreement with the ob¬ 
servations of R 136 and NGC 3125-A1. We note that the intrinsic 
continuum flux estimated by lYacca et afl d 19951) is even lower 
(log(Luv) ~ 37.33), which may be attributed to uncertainties in 
the UV extinction and the precise pointing of their observations. 

We conclude that the Hen emission of very massive WNh 
stars at Z » Zlmc is probably strong enough to explain the ob¬ 
served strength of the He n emission of NGC 3125-A1 and R136. 
In the following, we use the approximate average luminosity of 
the seven WNh stars in R 136 to define a representative luminos¬ 
ity of log(LwN5h/To) = 6.7. With the numbers above, we then 
obtain 


log(L H en/(erg/s)) = log(A W N5h) + 37.7, (1) 

where iVwNSh denotes the equivalent number of WN5h stars. 

The observed Hen line widths (FW HM) lie in the range 
of 1000-1400km/s for NGC 3125-A1 dChandar etal.l [2004; 
IWofford et alJ 12014) and ~ 2500 km/s for R136. The individ¬ 
ual WN5h stars in R136 have a similar spectral morphology 
as VFTS 682, for which we estimate FWHM ~ 2000km/s. The 
variations in the observed FWHM may indicate that the real stel¬ 
lar populations are more heterogeneous than we assumed in our 
simplified approach. 



1200 1300 1400 1500 1600 1700 1800 

X/K 

Fig. 1 . Synthetic UV spectra for different metallicities Z from wind 
models for WNh stars from lGrafener & Hamannl d 2008l) . The presented 
models are computed for a luminosity of log(L/L 0 ) = 6.3, a stel¬ 
lar temperature T* = 44.7 kK, and have very similar mass-loss rates 
(~ 1.8 x 1(T 5 M 0 yr~'). The given Hen line luminosities have only been 
determined for the narrow-emission components of the He ii line pro¬ 
files. Measurement errors for the given L HcI ,, FWHM, and EW probably 
amount to ~ ±0.05 dex (or ±10%). 


The He n line profiles in Fig.Qjconsist of a narrow emission 
component, a comparatively weak blue-shifted P-Cygni absorp¬ 
tion, and broad line wings that are caused by multiple electron 
scattering of line photons in the dense ionised wind. In addi¬ 
tion to the iron forest between 1300 and 1500 A, which is most 
prominent at high Z, the spectra show C iv, N iv, and N v features 
whose strength decreases for lower Z. We note that the mod¬ 
els have been computed for an N-enriched surface composition, 
which will affect the relative strength of the C iv vs. N iv and N v 
features. Furthermore, the resulting He n line fluxes will depend 
on the adopted stellar temperature T*. 

For low metallicities Z = 0.1 and 0.01, we find He n line lu¬ 
minosities of the order of loglLeen) ~ 37.5 for a given stellar 
luminosity of log (L/L 0 = 6.3). Adopting the same representa¬ 
tive stellar luminosity as in Eq.Q] we obtain 

log(L H en/(erg/s)) = logW W Nh,iowz) + 37.9, (2) 


3.2. He ii emission at iowZ 

To estimate L He n for WNh stars at low Z, we computed syn- 
thetic UV spectra using the wind models of lGrafener & Hamannl 
d2008l) . To obtain a more realistic FWHM, we used a line 
broadening velocity of 10 km/s instead of the artificially en¬ 
hanced value of lOOkm/s in the original work. In Fig. Q] we 
show three models with Z - Z 0 , 0.1 Z 0 , and 0 .01 Z 0 (cf. model 
grid 2 in Table 2 of lGrafener & Hamannll200'8h . For these mod¬ 
els, F e has been adjusted to reproduce very similar mass-loss 
rates for all three metallicitieQ As discussed in Sect. 12.21 the 
models with low Z have lower terminal wind velocities Vco and 
thus higher wind densities. This leads to stronger and narrower 
Hen emission at lower Z (cf. Fig.Q]). The resulting line fluxes 
and FWHM are indicated in Fig.Q] We note that the wind mod- 
els seem to underestimate v x at Z = Z Q (cf. the discussion in 
iGrafener & Hamamill2008l) . The important point is thus the pre¬ 
dicted qualitative behaviour, that is, the strong decrease of vs. 
Z. 

1 iDoran et~aD used the same lKroupal 120011 ) IMF as in this work. 

2 As the luminosities and surface abundances are kept fixed, this 
change in T e translates into different stellar masses M for each object. 


that is, a slightly higher value than the one that we obtained for 
LMC metallicity. 

While Lee u increases for lower Z, the strength of the ob¬ 
servable metal lines decreases. Most notably, all metal lines are 
substantially weaker than Hen and will not be detectable in 
the integrated spectr a of He n emitting galaxies, as observed by 
ICassata etaT. ' ( 201 3). 

We conclude that WNh stars at low Z can produce narrow 
He ii emission with similar or slightly stronger line fluxes than 
those discussed in Sect. l3.ll if the high required Eddington fac¬ 
tors r e can be maintained. Theoretically, high effective I c are 
indee d expected due to increa sed stellar rotation rates at low Z 
(e.g.. iMevnet & Maedeil[2()()2 ). Adopting similar relative num¬ 
bers and luminosities of WNh stars as in Sect. 13. II we would 
thus expect similar or slightly higher line fluxes and EWs than 
those discussed in Sect. l3.1l for young starbursts at low Z. 

We also note that rotating massive stars in extremely 
metal-poor environments (Z ~ 10 s ... 10 8 Z 0 ) are expected 
to reveal primary nitrogen with mass f r actions of about 1% 
at their surface (Yoon & Langerl 120051: IMevnet et alJ [2006). 
IGrafener & Hamami (2008) showed that such self-enriched ob¬ 
jects can develop very similar winds as WNh stars at low metal- 
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licity (roughly corresponding to Z ~ 0.02 Z Q ). We thus expect 
that even very early stellar generations can contain WNh stars 
with similar spectral properties_as di scus sed here. For Pop III 
stars with zero metal content. lEkstrom et ali ( 2008 ) showed that 
qualitatively different internal mixing properties can prevent 
strong self-enrichment, that is, their mass-loss propertie s may be 
different (but cf. also the discussion in lVink & de Koteril2005l) . 


4. Narrow Hen emission in star-forming galaxies at 
moderate redshifts 

A systematic study of star-forming galaxies with Heji emis¬ 
sion at moderate redshifts has been performed by ICassata et al.l 
(1201 3h . based on data from the VIMOS VLT Deep Survey 
(VVDS; iLe Fevre et aD l2005l). In the total VVDS sample of 
~ 45000 galaxy snectra JCassata et alJ found 277 star-forming 
galaxies with secure redshifts in the range 2 < z < 4.6, of which 
39 show He n A 1640 in emission. They divided this sample into 
four groups of 11 objects with narrow He n lines, 13 objects with 
broad Hen lines, 12 possible Hen emitters, and 3 AGN, to dis¬ 
tinguish between high-velocity outflows of WR stars, AGN, or 
supernova-driven winds, and the narrow nebular emission that is 
expected for Pop III stars. 

ICassata et al.l found that narrow and broad emitters show 
qualitatively different spectra. In particular, they found that only 
broad He n emitters (with FWHMfHe n) > 663 km/s) show in¬ 
dications of underlying P-Cygni type Civ 41550 emission. As 
we discussed in Sect. 12.21 this emission is most likely a signa¬ 
ture of stellar-wind emission from classical WR stars, predomi¬ 
nantly originating from WC subtypes at low Z. In the group of 
11 narrow emitters, they only found two such cases, in line with 
our prediction of very weak C iv emission for WNh stars at low 
Z. Our predicted FWHM(He ii) for such objects in Sect. 13. 21 is 
clearly within the limit for narrow emitters from ICassata et all 
even if we take possible uncertainties (up to a factor of ~ 2) into 
account. 

Excluding AGN, ICassata et al.l found He n line luminosities 
in the range log(L H en) ~ 40... 41.5, and EWs of « 1... 7A. They 
noted that, based on the predicted nebular Hen emission for 
Pop III objects from lSchaererl d2003l) , the star-formation rates re¬ 
quired to produce the observed line fluxes are about two orders 
of magnitude lower than those determined from the spectral en¬ 
ergy distribution (SFR se d). 

Here we furthermore note that the majority of narrow He u 
emitters is found at the bottom of the Hen luminosity distri¬ 
bution in a_very narrow range of log(Lneu) ~ 40...40.3 (cf. 
Fig. 10 in ICassata et alJ ( 2013 )). This low-luminosity group is 
additionally restricted to the lowest redshifts of z ~ 2.0 ...2.5. 
This implies that these objects are near the detection limit. Only 
three objects are located outside the low-luminosity group, with 
log(L H eu) ~ 41.0 ...41.5, and z = 2.5 ...4.0. 

The observed line fluxes for the low-luminosity group are 
remarkably similar to our expectations for massive young super¬ 
clusters with M c i ~ 10 6 M o at low Z. Clusters in this mass 
range are near the top end of the cluster initial mas s function 
(~ 2 x 10 6 Mm), as determined for instance hv iBastian l (l2008h for 
interacting galaxies and luminous IR galaxies. They may thus be 
the most luminous He n sources of this kind. For comparison, the 
observed continuum fluxes log(Luv) = log(LHen/EW(Hen)) ~ 
39.2 ...40.3 are similar to or higher than the Starburst99 pre¬ 
diction of log(Ljjv) ~ 39.04 for a cluster with 10 6 M Q and 
an age of 2-4 Myr at low Z = 0.0004. Test computations for 
galaxies with continuous SFRs at the same metallicity suggest 


Tuv £ 40.2 + log(SFR/(M Q yr 1 )), meaning that the host galax¬ 
ies could provide a substantial contribution to the observed con¬ 
tinuum. 

In our scenario, the continuum of a star-forming host galaxy 
could easily become so bright that it would prevent the detec¬ 
tion of He ii. The target selection by ICassata et al.l may thus be 
responsible for picking the most extreme and recent starburst 
events at low Z that are located in host galaxies with continua 
weak enough to enable the detection of He ii. The latter point 
places strong constraints on the continuum flux, which may be 
the reason for the narrow range of observed i-band magnitudes 
(23.45... 24.63 mag,) compared to the wide range of derived 
SFR se d ~ 1... 10 3 M q yr~' of Cassata et alJ (cf. their Table 2). 

5. Conclusions 

We discussed a new scenario for the origin of narrow He n emis¬ 
sion in star-forming galaxies at moderate redshifts, based on 
predictions of strong WR-type stellar winds with low terminal 
wind speeds for very massive stars (VMS) at low metallicities 
(Z ~ 0.01 Z 0 ). We estimated the expected Hen emission and 
found that it is in line with recent observations by ICassata et al.l 
(120131) . In our scenario the observed emission originates from a 
population of VMS in one or more young super-clusters located 
within a variety of star-forming galaxies at low Z. This scenario 
is in accordance with the observed line fluxes and EWs, the large 
observed spread of star-formation rates, and the different Civ 
profiles of the narrow and broad emitters. 

If our interpretation's correct, the narrow and broad Hen 
emitters identified bv ICassata et alJ both originate from massive 
stellar populations at low Z in distinct evolutionary stages, and 
there is no need to invoke the existence of Pop III stars at moder¬ 
ate redshifts to explain the observed narrow He n emission. The 
fact that the emission of VMS is largely neglected in current pop¬ 
ulation synthesis models will generally affect the interpretation 
of the integrated spectra of young stellar populations (cf. also 
IWofford et alllmfl) . 

Finally, we note that the same argumentation will hold for 
rotating VMS at extremely low Z (~ 10~ 5 ... 10 8 Z 0 .), for which 
similarly strong winds are predicted as a result of the self¬ 
enrichment with primary nitrogen. The He ii emission from such 
very early stellar generations of VMS may be detectable in fu¬ 
ture studies of star-forming galaxies at high redshifts with the 
JWST. 
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